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Crystalline order in dense packings of long polymers with a definite handedness is 
difficult to reconcile with the tendency of these chiral objects to twist and braid about 
each other. If the chirality is weak, the state of lowest energy is a triangular lattice of 
rigid rods. When the chirality is strong, however, screw dislocations proliferate, leading 
to either a tilt grain boundary phase or a new "moire state" with twisted bond order. In 
the latter case, polymer trajectories in the plane perpendicular to their average direction 
are described by iterated moire maps of remarkable complexity, reminiscent of dynamical 
systems. 
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1. Introduction 

A notable feature of biological materials is the profusion of long polymer molecules 
with a definite handedness. DNA, polypeptides (such as poly-7-benzyl-glutamate) and 
polysaccharides (such as xanthan) can all be synthesized with a preferred chirality. Long 
polymers in dense solution often crystallize into a hexagonal columnar phase, i.e., a lattice 
of rod-like objects with the cross section a triangular lattice. When the polymers are chiral 
this close packing into a triangular lattice competes with the tendency for the polymers 
to twist macroscopically [1] as in cholesteric liquid crystals. Similar to the twist grain 
boundary phase of chiral smectics [2], macroscopic chirality can proliferate when screw 
dislocations enter the crystal. Though the screw dislocations cost a finite energy, when the 
energy increase is offset by the decrease in free energy from twisting, the screw dislocations 
will penetrate the sample. If the chirality is weak, a defect free hexagonal columnar phase 
persists. 

In this note we discuss the effect of chirality on the hexagonal columnar phases [3] of 
long polymers [4] . We neglect for simplicity heterogeneity along the polymer backbones and 
work with a two component displacement field perpendicular to the local polymer direction. 
The usual chiral term relevant for cholesteric liquid crystals produces a polymer tilt grain 
boundary phase, similar to the smectic-^* phase [2] . We find, as well, an additional term 
in the free energy which favors the rotation of the bond order along the average polymer 
direction. This term leads to braided polymers with twisting describable by a sequence of 
moire patterns. The polymer braids can be described by a map from one perfect crystalline 
lattice to another perfect lattice rotated by a relative angle (j) and separated by honeycomb 
network of defects. When takes on certain lock-in angles, the free energy of the crystal 
distortion will go through a minimum. The resulting map creates moire patterns when the 
monomer locations are projected down the polymer axes. Upon iteration, corresponding 
to many periodically spaced defect walls, these moire maps exhibit complex, self-similar 
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behavior, reminiscent of chaotic dynamical systems. Little appears to be known about the 
properties of these iterated moire maps. 

In the next section we review the continuum elastic theory which applies to a polymer 
crystal. We show how the energy of dislocation defects in the crystal determines the critical 
values of the chiral coupling constants at which chirality first enters. We describe the 
smectic twist grain boundary state as well as its polymer analogue, the tilt grain boundary 
state. In the last section we describe a microscopic model of the moire phases, and discuss 
the lock-in angles which will produce moire maps and braided polymer crystals. Finally, 
in the appendix we prove that any rotation angle of a crystal which produces a non-trivial 
coincidence lattice (as do the special lock-in angles discussed above) is an irrational fraction 
of 27r. 



2. Statistical Mechanical Model and the Tilt Grain Boundary Phase 

In order to calculate averages of physical quantities, we integrate over the ensemble 
of all states, with each state weighted by the Boltzmann weight 

P[u] = Z-^ e^p{-F[u]/ksT} (2.1) 

where u{xi,X2,X3) is the two-dimensional displacement field of the polymers from a per- 
fect triangular lattice of rods, and Z is the partition function. The elastic free energy, 
appropriate for a crystalline state of polymers lying on average along the 2-axis, is 

F[u] = J ^ ^ Ksid^^Uif - 7V, xSfi- ^'d^e^^ (2.2) 

where Uij = ^{diUj + djUi), Sui = d^Ui, and 9q = \e^^i,pnf^di,Up = ^eijdiUj in the nematic 
(no = z) ground state. These fields measure respectively the local strain, tilt and bond 
order of the crystal. Here and throughout Roman indices run from 1 to 2 while Greek 
indices run from 1 to 3. Terms which are not invariant under spatial reflection are chiral 
and are not allowed in systems with mirror symmetry. The last two terms in (2.2) are the 
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only allowed couplings in a chiral polymer system to lowest order in the gradients of u. 
The term proportional to V_l x 5n favors a cholesteric-like twisting [3] of the local polymer 
direction along a pitch axis, perpendicular to the planes in which the polymers, on average, 
lie. The term proportional to 8^9^ favors a twisting of the bond order along the polymer 
direction. This term, and its effect on the polymer trajectories, can be described by moire 
maps. In the columnar crystal the two chiral terms are the same if d^di — did^. However, 
in the presence of crystal dislocations derivatives do not commute, V_l x5n = eijdidzUj ^ 
dzSijOiUj = SzOq. a dislocation line is characterized by the value of the line integral around 
this defect [5] 



where h is the Burgers vector. The Burgers vector must be an integer linear combination 
of lattice vectors. It points parallel to the defect axis for the screw dislocation shown in 
Figure 1. 

Various arrangements of dislocations lead to the twisting and braiding of polymers 
favored by 7 and 7'. Burgers vectors of dislocations in a hexagonal columnar phase lie in the 
x|/-plane, and there arc three generic types [3]: a screw dislocation, an edge dislocation with 
tangent along and an edge dislocation lying in the xy-plane. The latter defect requires 
aligned polymer ends which we can neglect if the polymers are very long. The remaining 
dislocations must lie in a plane spanned by their Burgers vector h and 2, which amounts 
to choosing dislocation complexions with axy = ocyxi where the dislocation density tensor 
a^i is the density of dislocations with tangents along the 7-direction with Burgers vectors 
pointing in the i direction [6] . 

Proceeding as in [5,6] we introduce a new field w^i = d^Ui away from any dislocations. 
The non-commutivity of the derivatives of u is represented by the dislocation density, 
^nu-ydi^Wji = —a^i. One can solve for w.yi in terms of the dislocation density a^i and find 
the equilibrium displacement field in the presence of crystal defects [4]. In terms of the 
nematic and bond order field non-commutivity of derivatives means 2dz9Q — Vj_ x 5n = 




(2.3) 
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— Trfo;]. The energy per unit length /g of a screw dislocation is finite, while that of an 
edge dislocation lying along the z direction diverges logarithmically with system size [7]. 
A screw dislocation is illustrated in Figure 1. Depending on how the dislocations are 
combined together a screw dislocation array can either produce a non-vanishing Vj_ x bn 
or dzQ&. Our predictions for the phase diagram as a function of 7 and 7' are summarized 
in Figure 6. 

If the chirality is strong and 7 7' we expect the polymer analogue of the Renn- 
Lubensky twist-grain-boundary state [2] of chiral smectic liquid crystals. Smectic liquid 
crystals are, in some sense, dual to the hexagonal columnar phases considered here [3]. 
Smectics are composed of sheets of rod-like molecules; they are crystalline in one direction 
(the direction perpendicular to the layers) instead of two. If the molecules are chiral these 
layers want to twist, thus disrupting the crystalline order. At low chiralities, macroscopic 
twist is suppressed. Above a certain threshold chirality, screw dislocations enter the sample, 
causing the smectic slabs to undergo, on average, a uniform rotation along a direction 
perpendicular to the layer normals. In Figure 2 we show the structure of this state. Polymer 
crystals, as well, can have a state characterized by a constant overall twist perpendicular to 
the polymer axes. The state consists of a parallel array of tilt grain boundaries, separated 
by cZ' along the rotation axis, bounding regions of perfect polymer crystals similar to the 
smectic planes in Figure 2. Each of these tilt grain boundaries (TGB) is composed of 
a parallel array of screw dislocations lying, say, in the x^-plane, pointing along the x- 
axis and uniformly spaced along z with spacing d. As illustrated in the upper inset of 
Figure 6 this dislocation texture causes a discrete rotation (/> = tan~-'^(6/(i) (with h = oq, 
corresponding to the shortest allowed Burgers vector) in the average polymer direction 
across the boundary. The spatial integral of Vj_x5n is non-zero, while the integral of dzOe 
vanishes. The TGB state appears when the chiral coupling 7 exceeds the critical value 
lc = fs/b [4]. 
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3. The Moire State 

The second chiral couphng 7' has no analogue in chiral smectics. To find a configura- 
tion of screw dislocations which exploits this form of chirality we search for a dislocation 
texture which produces a displacement depending only on z. The only texture which does 
not produce divergent elastic energy is one in which the only non-zero components of a^i 
are axx = otyy [6,4]. A honeycomb array of screw dislocations, on average, produces this 
dislocation texture. We now find that for a single such grain boundary x vanishes 
far from the boundary while dzQ^ does not. This sort of grain boundary thus causes a net 
twist of the hexatic bond order parameter while imposing no net VxX5n. When 7' » 7 
and the chirality is large, screw dislocations penetrate for 7' > 7c = 2/s/6. By choosing 
the rotation angle across the honeycomb dislocation network to produce a high density of 
coincidence lattice sites [8] we produce especially low strain energies across the boundary. 
The superposition of triangular polymer lattices below and above the boundary forms a 
moire pattern. The superposition of many such "moire sheets" along the 2-axis braids the 
polymers with deep minima in the energy at certain lock- in angles. Figure 3 illustrates 
the mapping of polymers across the moire plane. Polymers in the lower half-space (circles) 
must be connected to the closest available polymer in the upper half-space (crosses) to 
minimize bending energy. Note that the map has a discrete translational symmetry, in 
the sense that any coincidence site could be a center of rotation. Especially simple moire 
maps arise for the rotation angles 

0„ = 2 tan 

n = 1,2,.... Figure 3 shows these maps for n = 1, . . . , 4. It is shown in Appendix A 
that all such angles are irrational fractions of 27r [4] so that the structure never repeats 
upon iteration. Around each coincidence point there are n concentric rings of helical 
polymers. The lattice of coincidence points is also a triangular lattice, but with a spacing 
v/l + 3(2n+ 1)2/2, where ao is the original lattice constant. The geometrical 
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3(2n+ 1) 



(3.1) 



origin of such energetically preferred lock-in angles has no analogue in chiral smectics. 
The exact choice of lock-in angles and spacing between moire planes must be settled by 
detailed energetic calculations. 

Upon two iterations of the moire map separating three regions of polymer crystal, the 
first coincidence lattice is rotated with respect to the second coincidence lattice by precisely 
the angle of rotation 0^. Thus the composite coincidence lattice is the "coincidence lattice 
of coincidence lattices", with lattice constant a^/ao- Moire maps iterated p times lead to 
triangular composite coincidence lattices with spacing an{an/aoy~^ , i.e. to ever sparser 
lattices of fixed points with intricate fractal structure in between them. Figure 4 shows the 
projected polymer paths for a lock-in angle of a square lattice (with 0i = tan~-'^(3/4) 
36.9°) iterated p — 1, . . . , 4 times. For a square lattice an especially simple sequence of 
lock-in angles is given by 



-1 



2n + 1 



(3.2) 



_2n{n + l) 

leading to coincidence lattice spacings a^^^^^ = a^^n^ + (n -|- 1)^. 

In contrast to the TGB state polymers for p » 1 are highly entangled and wander far 
from straight line trajectories. The polymer configurations and the dislocations leading 
to them are shown in Figure 5 for the n = 1 moire map of a triangular lattice iterated 
nine times. Near the moire planes both $q and V_l x Sn are nonzero. The center polymer 
is a fixed point of all the maps. Any such fixed polymer has a halo of others twisting 
around it. In this special tube, the nematic order parameter takes on the texture of a 
double twist cylinder as found in the low-chirality limit of blue phases of cholesteric liquid 
crystals [9] The moire state takes advantage of both double twist energies and the new 
chiral coupling 7'. 

The phase diagram in Figure 6 summarizes our conclusions. As for chiral smectics [2] , 
it is possible to cast the theory in a form similar to the Ginzburg-Landau theory of Type 
II superconductors in a magnetic field [4]. The couplings 7 and 7' represent two distinct 
"magnetic fields" in this analogy. The hexagonal columnar phase expels all macroscopic 
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chirality. The two chiral couplings 7 and 7' cause screw dislocations to penetrate the crystal 
above critical strengths 7c and 7^, similar to the penetration of Abrikosov vortices above 
the lower critical field Hd- The TGB phase predicted here for chiral polymers is similar 
to that already observed experimentally in chiral smectics [2]. The braided moire state 
is qualitatively new, and its experimental observation would be of considerable interest. 
In Figure 7 we show 40 random polymers subject to 99 iterations of the n = 1 map for 
the triangular lattice, all near a central, fixed coincidence lattice point. Surprising little is 
known about the intricate trajectories produced by such iterated moire maps. These could 
be studied experimentally via neutron diffraction in hexagonal columnar crystals with a 
dilute concentration of deuterated polymer strands. Numerical studies of the local fractal 
dimension and Lyapunov exponent of these polymer trajectories are currently in progress. 
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Appendix A. Irrationality of the Lock-In Angles 

The special lock in angles that we have found for the triangular lattice are all irrational 
fractions of 27r. We will prove a more general theorem in this appendix that applies to any 
rotation which leads to a coincidence lattice with a finite lattice constant. 

Let R(^, A) = A' be the lattice generated by rotating a lattice of points A by ^ 
around a lattice site at the origin. Consider a regular lattice A 1, a set of ordered pairs 
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{{xi,yi), {x2, 1)2)1 ■ ■ ■}■ Suppose there is an angle 9* for which the intersection of Ai and 
R(^*, Ai) is a lattice (not equal, as in the pathological case Ai = R(^*, Ai), to the original 
lattice) with a finite lattice constant. These angles are special lock- in angles which keep 
a finite fraction of polymers from stretching. Let A" = R(n^*,Ai). Let the intersection 
of A5 and A} be the coincidence lattice A2, and similarly fl ^1 = ^T+i- Note that 
A^ = R(m6'*, A^). The intersection of A\ and A^ is also a coincidence lattice, but as 
Aj; and A^ are rotated by 9* from A^ and A}, respectively, their intersection, A2 will be 
rotated by 9* from A2. Thus, since the intersection of A?, A^ and Af is just the intersection 
of A2 and A2 the resulting coincidence lattice is the coincidence lattice of the coincidence 
lattices A2. 

In general, if we look at A^^.^ 



Thus A°_,_]^ is the coincidence lattice of two coincidence lattices A„. Therefore as we add 
layer after layer, consecutively rotating by 9* we will produce ever sparser lattices with 
growing lattice constants. 

Now assume that 9* is a rational fraction of 27r, so 9* = 27r|. Then, after q iterations of 
the rotation the lattices would repeat, or, in other words, A^ = A]*. This would mean that 
^q+i = 0^=0 -^1 ~ n^=o -^1 = since Af has already been included in the intersection of 
the first q lattices. This contradicts the conclusion of the preceding paragraph. Therefore 
it is impossible to have a rotation angle 9* which is both a rational fraction of 27r and 
produces a coincidence lattice with a non-trivial but finite lattice constant. Thus all lock- 
in angles must be irrational fractions of 27r, which is what we were to prove. 




(A.l) 
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Figure Captions 

Fig. 1. A single screw dislocation in a polymer crystal. The dark horizontal line is the 
screw dislocation. 

Fig. 2. The smectic twist grain boundary state. Each plane is a smectic layer and each 
stack is separated by a smectic Twist Grain Boundary. Figure provided courtesy 
of Tom Lubensky. 

Fig. 3. A single moire map for triangular lattices with n = 1,2, 3,4 (^^ = 21.8°, 13.2°, 
9.4°, 7.3° respectively). The crosses can be thought of as the ends of polymers 
in a perfect crystalline region, as can the circles. They join by connecting to the 
nearest polymer in order to reduce the elastic energy. The shaded lines are the 
screw dislocations which make up a honeycomb network. 

Fig. 4. The projected top view of a moire map on a square lattice with rotation angle 
tan~^(3/4). The four boxes show the projected polymer paths after p= 1, 2, 3, 4 
iterations. 

Fig. 5. The moire state. The thick tubes running in the z direction are polymers, while 
the dark lines are stacked honeycomb arrays of screw dislocations. The intersec- 
tion of these polymers with any constant z cross section away from the hexagonal 
defect arrays forms a perfect triangular lattice. 

Fig. 6. Phase diagram of a chiral polymer crystal. Insets are representative tilt (TGB) 
and moire grain boundaries. Shaded lines are screw dislocations. 

Fig. 7. A projected top view of 40 random polymers paths resulting from the moire map 
with n = 1 iterated 99 times. There is an exceptional fixed point of all 99 maps 
at the center. 
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